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Abstract 

We show that every weighted connected graph G contains as a subgraph a spanning tree 
into which the edges of G can be embedded with average stretch 0(log^ rtloglogn). Moreover, 
we show that this tree can be constructed in time 0(mlogri + nlog^ n) in general, and in time 
O(TOlogn) if the input graph is unweighted. The main ingredient in our construction is a novel 
graph decomposition technique. 

Our new algorithm can be immediately used to improve the running time of the recent solver 
for symmetric diagonally dominant linear systems of Spielman and Tcng from 

m2(°(^'°s"'°si°s")) to mlog^^^^n, 

and to 0(n log^ n log log n) when the system is planar. Our result can also be used to improve 
several earlier approximation algorithms that use low-stretch spanning trees. 

Categories and Subject Descriptors: F.2 [Theory of Computation]: Analysis of Algorithms 
and Problem Complexity 

General Terms: Algorithms, Theory. 

Keywords: Low-distortion embeddings, probabilistic tree metrics, low-stretch spanning trees. 



1 Introduction 

Let G = (y, E, w) be a weighted connected graph, where is a function from E into the positive 
reals. We define the length of each edge e G ii^ to be the reciprocal of its weight: 

d{e) = l/w{e). 
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(URI) administered by the Office of Naval Research under Grant N00014-01- 1-0795. The work was also partially 
supported by the Lynn and William Frankel Center for Computer Sciences. 

^Partially supported by NSF grant CCR-0324914. Part of this work was done at Yale University. 

*Partially supported by NSF grants CCR-0311430 and ITR CCR-0325630. 
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Given a spanning tree T of V, we define the distance in T between a pair of vertices u,v £ V, 
distT{u,v), to be the sum of the lengths of the edges on the unique path in T between u and v. 
We can then define the stretch^ of an edge {u, v) £ E to be 



stretchy (u, v) 



distT'(n, v) 



d{u, v) 



and the average stretch over all edges of E to be 



ave-stretchT(-E') = 7^ 2, stTetch.T{u,v). 
E — ^ 

{u,v)&E 



Alon, Karp, Peleg and West fl] proved that every weighted connected graph G = iV, E, w) of 
n vertices and m edges contains a spanning tree T such that 



and that there exists a collection r = {Ti, . . . , T^} of spanning trees of G and a probability distri- 
bution n over r such that for every edge e G E, 



The class of weighted graphs considered in this paper includes multi-graphs that may contain 
weighted self-loops and multiple weighted-edges between a pair of vertices. The consideration of 
multi-graphs is essential for several results (including some in jpy). 

The result of triggered the study of low-distortion embeddings into probabilistic tree metrics. 
Most notable in this context is the work of Bartal |51 1^1 which shows that if the requirement that 
the trees T be subgraphs of G is abandoned, then the upper bound of 1 can be improved by 
finding a tree whose distances approximate those in the original graph with average distortion 
0(logn • log log n). On the negative side, a lower bound of O(logn) is known for both scenarios 
01^. The gap left by Bartal was recently closed by Fakcharoenphol, Rao, and Talwar who 
have shown a tight upper bound of O(logn). 

However, some applications of graph-metric-approximation require trees that are subgraphs. 
Until now, no progress had been made on reducing the gap between the upper and lower bounds 
proved in 1 on the average stretch of subgraph spanning trees. The bound achieved in for general 
weighted graphs had been the best bound known for unweighted graphs, even for unweighted planar 
graphs. 

In this paper^, we significantly narrow this gap by improving the upper bound of from 
exp (O ( "v/log n log log n) ) to O(log^nloglogn). Specifically, we give an algorithm that for ev- 
ery weighted connected graph G = {V,E,w), constructs a spanning tree TOE that satisfies 
ave-stretchr = 0(log^ nloglogn). The running time of our algorithm is 0(m log n + nlog^n) 
for weighted graphs, and 0(m log n) for unweighted. Note that the input graph need not be simple 

^Our definition of the stretch diflers slightly from that used in jT|: distT{u,v)/dista{u,v), where distG{y,v) is the 
length of the shortest-path between u and v. See Subsection 11.11 for a discussion of the difference. 

^In the submitted version of this paper, we proved the weaker bound on average stretch of 0((log n log log n)^). 
The improvement in this paper comes from re-arranging the arithmetic in our analysis. Bartal _7_ has obtained a 
similar improvement by other means. 



ave- stretchT(£;) = exp (0(^1 og n log log n) 
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and its number of edges m can be much larger than (2). However, as proved in pP, it is enough to 
consider graphs with at most n(n + 1) edges. 

We begin by presenting a simpler algorithm that guarantees a weaker bound, ave-stretch7^(i?) = 
O(log^n). As a consequence of the result in |lj that the existence of a spanning tree with aver- 
age stretch /(n) for every weighted graph implies the existence of a distribution of spanning trees 
in which every edge has expected stretch /(n), our result implies that for every weighted con- 
nected graph G = (V^E^w) there exists a probability distribution 11 over a set r = {Ti, . . . ,Th} 
of spanning trees {T ^ E for every T € t) such that for every e € E, 'Et^u [stretch7^(e)] = 
0(log^ nloglogn). Furthermore, our algorithm itself can be adapted to produce a probability dis- 
tribution n that guarantees a slightly weaker bound of 0(log^ n) in time 0{m ■ log^ n). So far, we 
have not yet been able to verify whether our algorithm can be adapted to produce the bound of 
0(log^ n ■ log log n) within similar time limits. 

1.1 Applications 

For some of the applications listed below it is essential to define the stretch of an edge (n, v) £ E 
as in 1 , namely, stretchy (n, = distxiu, v)/dist g{u,v). The algorithms presented in this paper 
can be adapted to handle this alternative definition for stretch simply by assigning new weight 
w'{u,v) = l/distG{u,v) to every edge {u,v) G E (the lengths of the edges remain unchanged). 
Observe that the new weights can be computed in a preprocessing stage independently of the 
algorithms themselves, but the time required for this computation may dominant the running time 
of the algorithms. 

1.1.1 Solving Linear Systems 

Boman and Hendrickson 8 were the first to realize that low-stretch spanning trees could be used 
to solve symmetric diagonally dominant linear systems. They applied the spanning trees of |^ to 
design solvers that run in time 

^3/220(Vlognloglogn) log(l/e), 

where e is the precision of the solution. Spielman and Teng j21j improved their results to 

^20(Vlognloglogn) log(l/e). 

Unfortunately, the trees produced by the algorithms of Bartal E] and Fakcharoenphol, Rao, 
and Talwar cannot be used to improve these linear solvers, and it is currently not known 
whether it is possible to solve linear systems efficiently using trees that are not subgraphs. 

By applying the low-stretch spanning trees developed in this paper, we can reduce the time for 
solving these linear systems to 

mlog*^^^-* nlog(l/e), 

and to 0(n log^ n log log n log (1/e)) when the systems are planar. Applying a recent reduction of 
Boman, Hendrickson and Vavasis [H], one obtains a 0(nlog^ nloglognlog(l/e)) time algorithm 
for solving the linear systems that arise when applying the finite element method to solve two- 
dimensional elliptic partial differential equations. 
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1.1.2 Alon-Karp-Peleg-West Game 

Alon, Karp, Peleg and West constructed low-stretch spanning trees to upper-bound the value of 
a zero-sum two-player game that arose in their analysis of an algorithm for the /c-server problem: 
at each turn, the tree player chooses a spanning tree T and the edge player chooses an edge e & E, 
simultaneously. The payoff to the edge player is if e G T and stretchy (e) -|- 1 otherwise. They 
showed that if every n- vertex weighted connected graph G has a spanning tree T of average stretch 
/(n), then the value of this game is at most f{n) + 1. Our new result lowers the bound on the 
value of this graph-theoretical game from exp (0(\/log n log log n)) to O (log^ n log log n). 

1.1.3 MCT Approximation 

Our result can be used to improve drastically the upper bound on the approximability of the mini- 
mum communication cost spanning tree (henceforth, MCT) problem. This problem was introduced 
in [m, and is hsted as [ND7] in [H] and [TO]. 

The instance of this problem is a weighted graph G = {V, E, w), and a matrix {r(u, v) \ u,v £ V} 
of nonnegative requirements. The goal is to construct a spanning tree T that minimizes c(T) = 
T.u,vGV '^C^' ^) ■ distriu, v). 

Peleg and Reshef ^H] developed a 2*^(^'°s"''°s approximation algorithm for the MCT prob- 
lem on metrics using the result of A similar approximation ratio can be achieved for arbitrary 
graphs. Therefore our result can be used to produce an efficient 0(log^ nloglogn) approximation 
algorithm for the MCT problem on arbitrary graphs. 

1.1.4 Message- Passing Model 

Embeddings into probabilistic tree metrics have been extremely useful in the context of approxi- 
mation algorithms (to mention a few: buy-at-bulk network design [^j, graph Steiner problem jl3j . 
covering Steiner problem [121 )• However, it is not clear that these algorithms can be implemented 
in the message-passing model of distributed computing (see ^Hl)- In this model, every vertex of 
the input graph hosts a processor, and the processors communicate over the edges of the graph. 

Consequently, in this model executing an algorithm that starts by constructing a non-subgraph 
spanning tree of the network, and then solves a problem whose instance is this tree is very problem- 
atic, since direct communication over the links of this "virtual" tree is impossible. This difficulty 
disappears if the tree in this scheme is a subgraph of the graph. We believe that our result will 
enable the adaptation of the these approximation algorithms to the message-passing model. 

1.2 Our Techniques 

We build our low-stretch spanning trees by recursively applying a new graph decomposition that 
we call a star- decomposition. A star-decomposition of a graph is a partition of the vertices into sets 
that are connected into a star: a central set is connected to each other set by a single edge (see 
Figure^. We show how to find star-decompositions that do not cut too many short edges and 
such that the radius of the graph induced by the star decomposition is not much larger than the 
radius of the original graph. 

Our algorithm for finding a low-cost star-decomposition applies a generalization of the ball- 
growing technique of Awerbuch to grow cones, where the cone at a vertex x induced by a set of 
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vertices S is the set of vertices whose shortest path to S goes through x. 
1.3 The Structure of the Paper 

In Section|2 we define our notation. In Section|3J we introduce the star decomposition of a weighted 
connected graph. We then show how to use this decomposition to construct a subgraph spanning 
tree with average stretch O(log^n). In Section |1J we present our star decomposition algorithm. 
In Section we refine our construction and improve the average stretch to O (log^ n log log n) . 
Finally, we conclude the paper in Section 6 and list some open questions. 

2 Preliminaries 

Throughout the paper, we assume that the input graph is a weighted connected multi-graph G = 
{V,E,w), where w is a weight function from E to the positive reals. Unless stated otherwise, we 
let n and m denote the number of vertices and the number of edges in the graph, respectively. The 
length of an edge e £ E is defined as the reciprocal of its weight, denoted by d{e) = l/w^e). 

For two vertices u,v £ V, we define dist(ii,t;) to be the length of the shortest path between u and 
V in E. We write distG{u,v) to emphasize that the distance is in the graph G. 

For a set of vertices, 5 C 1/, G{S) is the subgraph induced by vertices in S. We write dist5(n, 
instead of distG'(5')(n, u ) when G is understood. E{S) is the set of edges with both endpoints in S. 

For S,T <^V, E{S, T) is the set of edges with one endpoint in S and the other in T. 

The boundary of a set S, denoted d (S), is the set of edges with exactly one endpoint in S. 

For a vertex x £ V, distG'(x, S) is the length of the shortest path between x and a vertex in S. 

A multiway partition of F is a collection of pairwise-disjoint sets {Vi, . . . , T4} such that |J^ Vi = V. 
The boundary of a multiway partition, denoted d (Vi, . . . , Vk), is the set of edges with endpoints in 
different sets in the partition. 

The volume of a set F of edges, denoted vol (F), is the size of the set \F\. 

The cost of a set F of edges, denoted cost (F), is the sum of the weights of the edges in F. 

The volume of a set S of vertices, denoted vol (S), is the number of edges with at least one endpoint 
in S. 

The ball of radius r around a vertex v, denoted B{r,v), is the set of vertices of distance at most r 
from V. 

The ball shell of radius r around a vertex v, denoted BS(r, u), is the set of vertices right outside 
B{r,v), that is, BS(r, u) consists of every vertex u £V — B{r,v) with a neighbor w £ B{r,v) such 
that dist(v, u) = dist(w, w) + d{u, w). 

For V £V , radc {v) is the smallest r such that every vertex of G is within distance r from v. For 
a set of vertices 5 C 1/, we write rads {v) instead of radc(5) {v) when G is understood. 
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3 Spanning- Trees of O(log^n) Stretch 



We present our first algorithm that generates a spanning tree with average stretch O(log'^n). We 
first state the properties of the graph decomposition algorithm at the heart of our construction. 
We then present the construction and analysis of the low-stretch spanning trees. We defer the 
description of the graph decomposition algorithm and its analysis to Section^ 

3.1 Low-Cost Star-Decomposition 

Definition 3.1 (Star-Decomposition). A multiway partition {Vq, . . . , V^} is a star- decomposition 
of a weighted connected graph G with center xq £ V (see Figure ^ if G Vq and 

1. for all < i < A;, the subgraph induced on Vi is connected, and 

2. for all i > 1, Vi contains an anchor vertex Xj that is connected to a vertex yj € Vq by an edge 
{xi,yi) G E. We call the edge {xi,yi) the bridge between Vq and Vi. 

Let r = radG(xo), and rj = rady^ (xi) for each < i < A;. For 5,e < 1/2, a star-decomposition 
{Vo, . . . , Vfc} is a {5, e)-star-decomposition if 

a. 5r < vq < (1 — 5)r, and 

b. ro + d{xi,yi) + < (1 + e)r, for each i > 1. 

The cost of the star-decomposition is cost {d (Vq, . . . , Vfe)). 




Figure 1: Star Decomposition. 

Note that if {Vq, . . . , Vk} is a {5, e) star-decomposition of G, then the graph consisting of the 
union of the induced subgraphs on Vo,...,Vk and the bridge edges {yk,Xk) has radius at most 
(1 + e) times the radius of the original graph. 

In Section m we present an algorithm StarDecomp that satisfies the following cost guarantee. 
Let X = (xi, ...,Xk) and y = {yi, . . .,yk)- 

Lemma 3.2 (Low-Cost Star Decomposition). Let G = {V, E, w) be a connected weighted graph 
and let xq be a vertex in V. Then for every positive e < 1/2, 

{{Vo, . . . ,Vk} ,x,y) = starDecomp(G, xq, 1/3, e) , 
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in time 0{m + n\ogn), returns a {1/3, e)- star- decomposition of G with center xq of cost 



cost{d{Vo,...,Vk)) < 



6mlog2(?7i + 1) 



e • radc (xq) 



On unweighted graphs, the running time is 0{m). 
3.2 A Divide-and-Conquer Algorithm 

The basic idea of our algorithm is to use low-cost star-decomposition in a divide-and-conquer 
(recursive) algorithm to construct a spanning tree. We use n (respectively, fh) rather than n (resp., 
m) to distinguish the number of vertices (resp., number of edges) in the original graph, input to 
the first recursive invocation, from that of the graph input to the current one. 

Before invoking our algorithm we apply a linear-time transformation from 1 that transforms 
the graph into one with at most n(n + l) edges (recall that a multi-graph of n vertices may have an 
arbitrary number of edges), and such that the average-stretch of the spanning tree on the original 
graph will be at most twice the average-stretch on this graph. 

We begin by showing how to construct low-stretch spanning trees in the case that all edges 
have length 1. In particular, we use the fact that in this case the cost of a set of edges equals the 
number of edges in the set. 

Fixa = (21og4/3(n + 6))-i. 

T = UnweightedLowStretchTree(G, xq). 

0. If |y| < 2, return G. (If G contains multiple edges, return a single copy.) 

1. Set p = radfj (a^o)- 

2. {{Vo, . . . , Vfc.} , X, y) = StarDecomp(G, xq, 1/3, a) 

3. For < i < k, 

set Ti = UiiweightedLowStretchTree(G(Vi), Xj) 



Theorem 3.3 (Unweighted). Let G = {V,E) be an unweighted connected graph and let xq be a 

vertex in V . Then 

T = UnweightedLowStretchTree(G', Xq) , 

in time 0{fhlogn) , returns a spanning tree of G satisfying 



Proof. For our analysis, we define a family of graphs that converges to T. For a graph G, we let 



4. SetT = \J^TiU\Ji{y^,x^). 



radT (xo) < \/e • radc (xq) 



(1) 



and 



ave- stretchy {E) < O(log^m) . 



(2) 



({Fo, . . . ,Vk} ,x,y) = StarDecomp(G, xq, 1/3, a) 
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and recursively define 



Ro{G) = G and Rt{G) = |J(yi, rr,) U J i?t_i(G(V-)) . 

i i 

The graph Rt{G) is what one would obtain if we forced UnweightedLowStretchTree to return its 
input graph after t levels of recursion. 

Because for all n > 0, (21og4/3(n + 6))"^ < 1/12, we have (2/3 + a) < 3/4. Thus, the depth 
of the recursion in UnweightedLowStretchTree is at most log4/3n, and we have R\og^f.^n{G) = T. 
One can prove by induction that, for every t > 0, 

rad^^(G) (xo) < (1 + a)*radG (xq) . 

The claim in (^Q) now follows from (1 + a)^°^4/3^ < To prove the claim in we note that 

stretchT(u, u) < (distT(2;o, + distT(a;o, f )) (3) 

{u,v)(id(Vo,...yk) {«,D)e9(Vb,...,Vfc) 

< ^ 2V^.radG(xo), by dU (4) 
{«,t,)e9(Vb,...,Vfc) 

r- /6mlog2(wi + 1)\ , , 

< 2ve • radc \xq) — - — r — , by Lemma I3T2I 

V a-radG(xo) / 

Applying this inequality to all graphs at all log4/3n levels of the recursion, we obtain 

stretchy (ti, u) < 24 -v/e mlog2 mlog4/3 n log4/3(n + 6) = O(mlog^m) . 

□ 

We now extend our algorithm and proof to general weighted connected graphs. We begin 
by pointing out a subtle difference between general and unit- weight graphs. In our analysis of 
UnweightedLowStretchTree, we used the facts that radc {xq) < n and that each edge length is 1 
to show that the depth of recursion is at most log4/3 general weighted graphs, the ratio of 

radc (xq) to the length of the shortest edge can be arbitrarily large. Thus, the recursion can be 
very deep. To compensate, we will contract all edges that are significantly shorter than the radius 
of their component. In this way, we will guarantee that each edge is only active in a logarithmic 
number of iterations. 

Let e = {u,v) be an edge in G = {V,E,w). The contraction of e results in a new graph by 
identifying u and f as a new vertex whose neighbors are the union of the neighbors of u and v. All 
self-loops created by the contraction are discarded. We refer to u and v as the preimage of the new 
vertex. 

We now state and analyze our algorithm for general weighted graphs. 
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Fix/3 = (21og4/3(n + 32))-i. 

T = LowStretchTree(G = {V, E, w),Xo). 

0. If < 2, return G. (If G contains multiple edges, return the shortest copy.) 

1. Set p = radc (xq)- 

2. Let G = {V, E) be the graph obtained by contracting all edges in G of length less than 
(5p/n. 

3. (^1^0, • • • , , 5, = StarDecomp(G, xq, 1/3, /?)• 

4. For each z, let Vi be the preimage under the contraction of step 2 of vertices in V^, and 
{xi,yi) G Vo X be the edge of shortest length for which Xi is a preimage of and yi 
is a preimage of y^. 

5. For < f < /c, set Tj = LowStretchTree(G(Vi), Xj) 

6. Setr = UT,uU(y„x,). 

In what follows, we refer to the graph G, obtained by contracting some of the edges of the graph 
G, as the edge- contracted graph. 

Theorem 3.4 (Low-Stretch Spanning Tree). Let G = (y,E,w) be a weighted connected graph 
and let xq be a vertex in V . Then 

r = LowStretchTree(G, Xo) , 

in time 0(mlogn + nlog^ n), returns a spanning tree of G satisfying 

radr (xo) < 2-v/e • radg (xq) (5) 

and 

ave-stretchT{E) = 0{log^ rh) (6) 

Proof. We first establish notation similar to that used in the proof of Theorem 13.31 Our first step 
is to define a procedure, SD, that captures the action of the algorithm in steps 2 through 4. We 
then define Ro{G) = G and 

Rt{G) = \J{yi,Xi)u\jRt-i{G{V,)), 

i i 

where 

({Vo, . . . , Vfc} ,xi, . . . ,Xfc,yi, . . . ,yj.) = SD(G,xo,l/3,/3). 

We now prove the claim in ©. Let p = radG'(2;o). Let t = 2 log4/3 (n + 32) Pt — 

rad^j(Q) (xq). Each contracted edge is of length at most f3p/n, and every path in the graph G{Vi) 
contains at most n contracted edges, hence the insertion of the contracted edges to G{Vi) increases 
its radius by an additive factor of at most Pp. Since (21og4/3(n + 32))-i < 1/24 for every n > 0, it 
follows that 2/3 + 2 /? < 3/4. Therefore, following the proof of Theorem 13. 3( we can show that pt 
is at most -^/e • radc (xq). 
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We know that each component of G that remains after t levels of the recursion has radius at 
most p(3/4)* < p/n^. We may also assume by induction that for the graph induced on each of 
these components, LowStretchTree outputs a tree of radius at most 2y/e{p/v?). As there are at 
most n of these components, we know that the tree returned by the algorithm has radius at most 

Vc/O + n X 2^(p/n^) < 2y/ep, 

for n > 2. 

We now turn to the claim in @, the bound on the stretch. In this part, we let Et<^E denote 
the set of edges that are present at recursion depth t. That is, their endpoints are not identified 
by the contraction of short edges in step 2, and their endpoints remain in the same component. 
We now observe that no edge can be present at more than log4/3 iC^^/P) + 1) recursion depths. 
To see this, consider an edge {u, v) and let t be the first recursion level for which the edge is in 
Ef. Let pt be the radius of the component in which the edge lies at that time. As u and v are not 
identified under contraction, they are at distance at least Ppt/n from each other. (This argument 
can be easily verified, although the condition for edge contraction depends on the length of the 
edge rather than on the distance between its endpoints.) If u and v are still in the same graph on 
recursion level t + log4/3 ((2^//^) + 1); then the radius of this graph is at most /?t/((2n//3) + 1), 
thus its diameter is strictly less than Ppt/n, in contradiction to the distance between u and v. 

Similarly to the way that Yl(u v)£d{Vo v^) ^) upper-bounded in inequalities ©- 

Q in the proof of Theorem 13.31 it follows that the total contribution to the stretch at depth t is 
at most 

0(vol (Et) log^ in). 

Thus, the sum of the stretches over all recursion depths is 

O (vol (Et) log^ m) = O [m log^ fh) . 

t 

We now analyze the running time of the algorithm. On each recursion level, the dominant cost 
is that of performing the StarDecomp operations on each edge-contracted graph. Let Vt denote 
the set of vertices in all edge-contracted graphs on recursion level t. Then the total cost of the 
StarDecomp operations on recursion level t is at most 0(|St| -|- | Vi| log | Vt|). We will prove soon 
that \ Vt\ = O(nlogn), and as |ii^t| = O(mlogm), it follows that the total running time is 
0(m log n -|- n log^ n). Note that for unweighted graphs G, the running time is only 0(m log n). □ 

The following lemma shows that even though the number of recursion levels can be very large, 
the overall number of vertices in edge-contracted graphs appearing on different recursion levels is 
at most 0(n log n). This lemma is used only for the analysis of the running time of our algorithm; 
a reader interested only in the existential bound can skip it. 

Lemma 3.5. Let Vt be the set of vertices appearing in edge- contracted graphs on recursion level t. 
Then Y.t\Vt\ = 0(n log n). 

Proof. Consider an edge-contracted graph G = {V , E) on recursion level t and let x be a vertex in 
V . The vertex x was formed as a result of a number of edge contractions (maybe 0). Consequently, 
X can be viewed as a set of all original vertices that were identified together to form it, i.e., x can be 
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viewed as a super-vertex. Let x{^) denote the set of original vertices that were identified together 
to form the super-vertex v £ V. 

We claim that for every super- vertex x £ Vt+i, there exist a super- vertex y £ Vt such that 
xi^) ^ x{y)- To prove it, note that every graph on recursion level t -|- 1 corresponds to a single 
component of a star decomposition on recursion level t. Moreover, an edge that was contracted 
on recursion level t + 1 must have been contracted on recursion level t as well. Therefore we can 
consider a directed forest J-' in which every node at depth t, corresponds to some super-vertex in 
Vt, and an edge leads from a node y at depth t to a node x at depth t + 1, if xi^) ^ xiv)- Note that 
the roots of J- correspond to super-vertices on recursion level and the leaves of correspond to 
the original vertices of the graph G. 

In the proof of Theorem 13. 4( we showed that every edge is present on O(logn) recursion levels. 
Following a similar line of arguments, one can show that every super- vertex is present on O(logn) 
(before it decomposes to smaller super- vertices, each contains a subset of its vertices). Since there 
are n vertices in the original graph G, there are n leaves in J^. Therefore the overall number of 
nodes in the directed forest is O(nlogn), and the bound on J2t 1^1 holds. 

□ 



4 Star decomposition 

Our star decomposition algorithm exploits two algorithms for growing sets. The first, BallCut, is 
the standard ball growing technique introduced by Awerbuch 4 , and was the basis of the algorithm 
of The second, ConeCut, is a generalization of ball growing to cones. So that we can analyze 
this second algorithm, we abstract the analysis of ball growing from the works of |16l Ol[T7|. Instead 
of nested balls, we consider concentric systems, which we now define. 

Definition 4.1 (Concentric System). A concentric system in a weighted graph G = {V,E,w) 
is a family of vertex sets C = {Lr QV:r£ M"*" U {0}} such that 

1. Lo + 0, 

2. Lr ^ Lr' for all r < r', and 

3. if a vertex u £ Lr and (n, v) is an edge in E then v £ Lr+d{u,v)- 

For example, for any vertex x £ V, the set of balls {B{r, x)} is a concentric system. The radius 
of a concentric system C is radius{C) = min{r : = V}. For each vertex v £V , we define ||f 
to be the smallest r such that v £ Lr. 

Lemma 4.2 (Concentric System Cutting). Let G = {V,E,w) be a connected weighted graph 
and let C = {Lr} be a concentric system. For every two reals < A < A', there exists a real 
r £ [A, A') such that 



vol (Lr) + T 

cost[o[Lr)) < — r max 

A — A 



l,log2 



m + r 



Y0l{E{Lx))+T 



where m = \E\ and 

^^fl ifyol{E{Lx)) = 
I otherwise. 
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Proof. Note that rescahng terms does not effect the statement of the lemma. For example, if all 
the weights are doubled, then the costs are doubled but the distances are halved. Moreover, we 
may assume that A' < radius{C), since otherwise, choosing r = radius{C) implies that d (Ly.) = 
and the claim holds trivially. 

Let Tj = \\vi\\, and assume that the vertices are ordered so that ri < r2 < • • • < r„. We may 
now assume without loss of generality that each edge in the graph has minimal length. That is, 
an edge from vertex i to vertex j has length | r j — rj\. The reason we may make this assumption 
is that it only increases the weights of edges, making our lemma strictly more difficult to prove. 
(Recall that the weight of an edge is the reciprocal of its length.) 

Let Bi = Lr^ ■ Our proof will make critical use of a quantity /ij , which is defined to be 



Tk — fj 
ivj,Vk)GE:j<i<k 

That is, fii sums the edges inside B^, proportionally counting edges that are split by the boundary 
of the ball. The two properties of /Xj that we exploit are 

m+i = + cost {d {Bi)) (rj+i - ri), (7) 

and 

T + vol{E{Bi)) < fii<T + vo\{Bi). (8) 
The equality follows from the definition by a straight-forward calculation, as 



= T + vo\{E{Bi+i)) -vol{{{vj,Vi+i) e E \ j <i}) + 

and 



rk 

ivj,v^)eE:j<i<k 



cost {d (Bi)) (ri+i - ri) = ^ 



n+i - Vi 



rk — 

{vj,Vk)eE:j<i<k ■' 

Choose a and b so that ra-i < A < and rf, < A' < Let = A' — A. We first consider 

the trivial case in which b < a. In that case, there is no vertex whose distance from vq is between 
A and A'. Thus every edge crossing i(A+A')/2 has length at least v, and therefore cost at most l/v. 
Therefore, by setting r = (A + A')/2, we obtain 

cost (d (Lr)) < vol {d (Lr)) ^ < vol (Lr) ^ , 

establishing the lemma in this case. 
We now define 

, f m + T 
ri = log2 



wo\{E{Ba-i))+T^ 

Note that Ba-i = Lx, by the choice of a. A similarly trivial case is when [a,b] is non-empty, and 
where there exists an i G [a — 1, 6] such that 



V 

r-i+i - r-j > - 
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In this case, every edge in d (Bi) has cost at most r]/u, and by choosing r to be max{ri,A}, we 
satisfy 

cost (d(Lr)) <\d(Lr)\- <Y0\(Lr)- , 

V V 

hence, the lemma is established in this case. 

In the remaining case that the set [a, 6] is non-empty and for alH G [a — 1, 6], 

V 

Ti+i - < - , (9) 

T] 

we will prove that there exists an z G [a — 1, 6] such that 

cost {d{Bi)) < ^ , 

hence, by choosing r = maxjrj. A}, the lemma is established due to ©. 
Assume by way of contradiction that 

cost id{Bi)) > fj,ir]/u 

for alH G [a - 1, 6]. It follows, by 0, that 

fXi+i > m+ /Uj(rj+i - ri)rj/i^ 

for all z G [a — 1, 6], which implies 

b 



i=a—l 
b 

> n 2(('^^+i , by (jH) and since 1 + x > 2"^ for every < x < 1 

i=a—l 

> fIa-l-{{m + T)/{Y0l{E{Ba-l))+T)) 

> m + T ,hy 

which is a contradiction. □ 

An analysis of the following standard ball growing algorithm follows immediately by applying 
Lemma 14.21 to the concentric system {B{r,x)}. 



r = BallCut(G, xo,p, S) 

1. Set r = 5 p. 

2. While cost {d {B{r, xq))) > log2(m + 1), 

a. Find the vertex v B{r,XQ) that minimizes dist(xo,w) and set r = dist(xo,'y). 
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Corollary 4.3 (Weighted Ball Cutting). Let G = {V,E,w) he a connected weighted graph, let 
X G V , p = radc (x), r = BallCut(G, xq, p, 1/3), and Vq = B{r, x). Then p/3 < r < 2 p/3 and 

^st om, < 3(v°iw) + i)i°g.(ii^i + i) , 

P 

We now examine the concentric system that enables us to construct Vi, . . . ,Vk in Lemma 13.21 

Definition 4.4 (Ideals and Cones). For any weighted graph G = {V,E,w) and S '^V, the set 
of forward edges induced by S is 

F{S) = {{u v) : (n, v) G E, dist(u, S) + d{u, v) = dist(t;, S)]. 

For a vertex v gV , the ideal of v induced by S, denoted Isiv), is the set of vertices reachable from 
V by directed edges in F{S), including v itself. 

For a vertex v £ V, the cone of width I around v induced by S, denoted Gs{l,v), is the set of 
vertices in V that can be reached from t> by a path, the sum of the lengths of whose edges e that 
do not belong to F{S) is at most /. Clearly, Cs{0,v) = Is{v) for all v £V. 

That is, Isiv) is the set of vertices that have shortest paths to S that intersect v. Also, 
u £ Gs{l,v) if there exist ao, . . . , Uk-i and bi,...,bk such that oq = v, bk = u, S /^(aj), 
{bi,ai) G E, and 

^d{bi,ai) < /. 

i 

We now establish that these cones form concentric systems. 

Proposition 4.5 (Cones are concentric). Let G = {V, E, w) be a weighted graph and let S C y. 
Then for all v GV, {Cs{l,v)}i^ is a concentric system in G. 

Proof. Clearly, Cs{l,v) C Cs{l',v) if I < I'. Moreover, suppose u £ Cs{l,v) and {u,w) £ E. Then 
if {u ^ w) £ F, then w £ Cs{l,v) as well. Otherwise, the path witnessing that u £ Cs{l,v) 
followed by the edge {u,w) to w is a witness that w £ Cs{l + d{u,w),v). □ 



r = ConeCut(G, v, A, A', S) 

1. Set r = A and if vol (^^(^^(A, v))) = 0, 

Set p = (vol (Cs(r, v)) + 1) log2(m + 1) 
otherwise. 

Set p = vol {Gs{r, v)) log2(m/vol {E{Cs{X, v))) 

2. While cost {d{Cs{r,v))) > p/{X' - A), 

a. Find the vertex w Gs{r,v) minimizing dist{w, Cs{r,v)) and set r 
dist{w, Gs{r,v)). 



r + 



Corollary 4.6 (Cone Cutting). Let G = (y,E,w) be a connected weighted graph, let v be a 
vertex in V and let S QV. Then for any two reals < A < A', ConeCut(G, f , A, A', S") returns a 
real r £ [A, A') such that 



^/f>fn f \\\ ^ vol (Cs(r, u)) +T 

cost[o [Cs[r,v))) < — max 

A — A 



l,log2 



m + T 



voliE{Gs{X,v)))+T 
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where m = \E\, and 

^ ^ r 1, ifyol{E{Cs{X,v)))=0 
\ 0, otherwise. 

We will use two other properties of the cones Cs{l, v): that we can bound their radius (Proposi- 
tion l4.7|) . and that their removal does not increase the radius of the resulting graph (Proposition l4.9j) . 

Proposition 4.7 (Radius of Cones). Consider a connected weighted graph G = {V,E,w), a 
vertex subset S and let = maxt,gy dist{v, S). Then, for every x £ S, 

radcs(;,x) (x) < + 21. 

Proof. Let u be a vertex in Cs{l, x), and let ao, • • • , Ok-i and 6i, . . . , 6^ be vertices such that oq = x, 
bk = u, bi+i £ Isiai), (bi,ai) £ E, and 

'^d{bi,ai) < I. 

i 

These vertices provide a path connecting x to u inside Cs{l,x) of length at most 

'^d{bi,ai) + ^dist(ai,6j+i). 

i i 

As the first term is at most /, we just need to bound the second term hy ^|J + 1. To do this, consider 
the distance of each of these vertices from S. We have the relations 

dist(6i+i, S*) = dist(ai, S) + dist(ai, bi+i) 
dist(ai, S) > dist(6j, S) — d{bi,ai), 

which imply that 

> dist{bk,S) >'Y^dist{ai,bi+i) - d{bi,ai) > dist(aj, - Z , 

as desired. □ 

In our proof, we actually use Proposition 14.81 which is a slight extension of Proposition 14.71 It's 
proof is similar. 

Proposition 4.8 (Radius of Cones, II). Consider a connected weighted graph G = {V,E,w), a 
vertex xq £ V and let p = radc (xq). Consider a real < p and let Vq = B{ro,xo), V = V — Vq 
and S = BS(ro,xo). Consider a vertex xi £ S and let ip = p — distG{xo,xi). Then the cones 
Cs{l,xi) in the graph G{V') satisfy 

Trades (i,xi) (xi) <4^ + 2l. 

Proposition 4.9 (Deleting Cones). Consider a connected weighted graph G = {V, E,w), a vertex 
subset S V , a vertex x £ S and a real I > and let V' = V — Gs{l,x), S' = S — Cs{l,x) and 
ijj = max^gy dist{v, S). Then 

max distv'iv, S') < ip . 



15 



Proof. Consider some v £ V . If the shortest path from v to S intersects Cs{l, x), then v G Cs{l, x). 
So, the shortest path from t; to S in y must he entirely in V'. □ 



The basic idea of StarDecomp is to first use BallCut to construct Vq and then repeatedly apply 
ConeCut to construct Vi, . . . ,Vk. 



{{Vo,...,Vk,x,y}) = StarDecomp(G = {V,E,w),xo,S,e) 

1. Set p = radfj (xq); Set tq = BallCut(G, xq, p, 6) and Vq = B{rQ, xq); 

2. Let S = BS(ro,a;o); 

3. Set G' = {V, E', w') = G{V - Vq), the weighted graph induced by V - Vq] 

4. Set {{Vi, Vk, x}) = ConeDecomp(G', 5, ep/2); 

5. For each i £ [1 : k], set yk to be a vertex in Vq such that {xk,yk) G E and is on a 
shortest path from xq to Xk- Set y = (yi, . . . , y^). 

{{Vi, . . . , Vfc, x}) = ConeDecomp(G, S, A) 

1. Set Go = G, Sq = S, and A; = 0. 

2. while Sk is not empty 

a. Set k = k + 1; Set x^ to be a vertex of S^-i; Set = ConeCut(Gfc_i, x^, 0, A, S^-i) 
h. Set Vk = Cs,_, irk,Xk); Set G^ = G(y - U^iVfe) and 5^ = - Vk- 

3. Set a; = (xi, . . . , x^). 



Proof of Lemma XH.iA Let p = radc {xq). By setting (5 = 1/3, Corollary 1^21 guarantees /o/3 < tq < 
(2/3)/9. Applying A = ep/2 and Propositions 14. 71 and 14.91 we can bound for every i, vq + d{xi,yi) + 
fi < p + 2A = p + ep. Thus StarDecomp(G, xq, 1/3, e) returns a (1/3, e)-star-decomposition with 
center xq. 

To bound the cost of the star-decomposition that the algorithm produces, we use Corollaries 
ESI and EH 

cost(8(y.))<^"+™'"^''»'°'^''"'+". and 

cost {e (v„ V - u^v;)) ^ni + ^oHVMn + i) 

for every 1 < i < /c, thus 

k 

cost (a (Vo , . . . , i^fc) ) < 5^ cost (y, , y - u^o^. 

< 2 1og.(m + l) j2i.oHV,) + l 

j=0 

6mlog2(m + 1) 



< 
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To implement StarDecomp in 0(m + n log n) time, we use a Fibonacci heap to implement steps 
(2) of BallCut and ConeCut. If the graph is unweighted, this can be replaced by a breadth- first 
search that requires 0{m) time. □ 



5 Improving the Stretch 

In this section, we improve the average stretch of the spanning tree to O (log^ n log log n) by intro- 
ducing a procedure ImpConeDecomp which refines ConeDecomp. This new cone decomposition trades 
off the volume of the cone against the cost of edges on its boundary (similar to Seymour ^20^). Our 
refined star decomposition algorithm ImpStarDecomp is identical to algorithm StarDecomp, except 
that it calls 

({^1, . . . , Vfc, x}) = ImpConeDecomp(G', 5, t, m) 
at Step 4, where t is a positive integer that will be defined soon. 



({Vi, . . . , Vfc, a;}) = ImpConeDecomp(G, S", A, t, fh) 

1. Set Go = G, Sq = S and j = 0. 

2. while Sj is not empty 



a. Set j = j + 1, set Xj to be a vertex of Sj^i, and set p = t — 1; 

b. while p > 

i. =ConeCut(G,_i,x,-, (i^^, (i^,5,„i); 

ii. if vol {E{Cs^_^{rj,Xj))) < Zt then exit the loop; else p = p - 1; 



2logP 

c. Set Vj = Csj^iirj, Xj) and set Gj = G{V — ul^-^Vi) and Sj = Sj-i — Vj] 
3. Set X = {xi, . . . , Xfc). 



Lemma 5.1 (Improved Low-Cost Star Decomp). Let G, xq and e be as in Lemma W^ t be 
a positive integer control parameter, and p = radc (xq). Then 

({Vo, . . . , Vfc}, X, y) = InipStarDecomp(G, xq, 1/3, e, t, fh) , 
in time 0{m + n\ogn), returns a {\/'i,e)- star- decomposition of G with center xq that satisfies 

cost{d{Vo)) < 6vol(V-o)log2(^ + l) 

p 

and for every index j G {1, 2, . . . , h} there exists p = p{j) € {0, 1, . . . , t — 1} such that 

(^(,, ,,\\ 4vol(y,)log(P+^)/*(m + l) , , 

cost [E [Vj , V - Uj^o^, j J < t ^ ^ ' (10) 

and unless p = 0, 

711 

vol(S(V,)) < J— . (11) 
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Proof. In what follows, we call p{j) the index-mapping of the vertex set Vj. We begin our proof 
by observing that < r j < ep/2 for every 1 < j < k. We can then show that {Vq, . . . , 14} is a 
(1/3, e)-star decomposition as we did in the proof of Lemma IX^ 

We now bound the cost of the decomposition. Clearly, the bound on cost(9(Vo)) remains 
unchanged from that proved in Lemma but here we bound vol (Vq) + 1 by 2vol (Vq). 

Below we will use A = ep/2 as specified in the algorithm. 

Fix an index j G {1, 2, . . . , k}, and let p = p[j) be the final value of variable p in the loop 
above (that is, the value of p when the execution left the loop while constructing Vj). Observe that 
p € {0, 1, . . . , f — 1}, and that unless the loop is aborted due top = Q, we have vol {EiVj)) < 
and inequality ((TT|) holds. 

For inequality ([T?H) . we split the discussion to two cases. First, consider the case p = t — 1. 
Then the inequality cost (^E (Vj,V - U^=o^i)) ^ (vol(V,) + l)log(m + l)(t/A) follows directly 
from Corollarv 14.61 and inequality l\W\i holds. 

Second, consider the case p < t — 1 and let r'^ be the value of the variable rj at the beginning of 
the last iteration of the loop (before the last invocation of Algorithm ConeCut). In this case, observe 

(as otherwise the 



'J^J — 2logP/* 7 



that at the beginning of the last iteration, vol \^E{Csj_j^{rj, Xj))j > ^jpg(p+i)/t ^ 
loop would have been aborted in the previous iteration). By Corollarv 14. H| 



costiE[V,,V-UU,V, < 



A/t 



X max 



l,log2 



m 



vol C5,_, 



(<-P-l)A ^. 



where Vj = Csj_i{rj,Xj). Since 



it follows that 



(t-p-2)A ^ ^, ^ (t-p-l)A 



^ol[E[Cs,_A^^^T^^^^)) > yol{E{Cs,_Ar'„x,))) > 



m 



Therefore 



and 



log(P+^)/*m > max 



l,log2 



m 




A/t 



□ 



Our improved algorithm ImpLowStretchTree(G, xq, t, ^ri), is identical to LowStretchTree ex- 
cept that in Step 3 it calls ImpStarDecomp(G, xq, 1/3, /3, t, m), and in Step 5 it calls 
ImpLowStretchTree(G(Vi), Xj, i, m). We set t = log log n throughout the execution of the algo- 
rithm. 
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Theorem 5.2 (Lower— Stretch Spanning Tree). Let G = {V,E,w) be a connected weighted 
graph and let xq be a vertex in V . Then 

T = InipLowStretchTree(G, xq, t, rn) , 

in time 0(mlogn + nlog^n), returns a spanning tree of G satisfying 

radr (xq) < 2i/e • radc (x) 

and 

ave- stretchy {E) = O (log^ n log log n) . 

Proof. The bound on the radius of the tree remains unchanged from that proved in Theorem 13.41 
We begin by defining a system of notations for the recursive process, assigning for every graph 
G = (y, E) input to some recursive invocation of Algorithm ImpLowStretchTree, a sequence (t{G) 
of non-negative integers. This is done as follows. If G is the original graph input to the first 
invocation of the recursive algorithm, than (t{G) is empty. Assuming that the halt condition of the 
recursion is not satisfied for G, the algorithm continues and some of the edges in E are contracted. 
Let G = (y, E) be the resulting graph. (Recall that we refer to G as the edge-contracted graph.) Let 
{Vb, 1^1, . . . , Vfc} be the star decomposition of G. Let V^- E ^ be the preimage under edge contraction 
oi Vj ^ V for every < j < A;. The graph GiVj) is assigned the sequence a{G{yj)) = (j{G) ■ j. 
Note that |o"(G)| = h implies that the graph G is input to the recursive algorithm on recursion level 
h. We warn the reader that the edge-contracted graph obtained from a graph assigned with the 
sequence a may have fewer edges than the edge-contracted graph obtained from the graph assigned 
with the sequence a-j, because the latter may contain edges that were contracted out in the former. 

We say that the edge e is present at recursion level /i if e is an edge in G which is the edge- 
contracted graph obtained from some graph G with |(t(G)| = h (that is, it was not contracted out). 
An edge e appears at the first level h at which it is present, and it disappears at the level at which it 
is present and its endpoints are separated by the star decomposition. If an edge appears at recursion 
level h, then a path connecting its endpoints was contracted on every recursion level smaller than 
h, and no such path will be contracted on any recursion level greater than h. Moreover, an edge is 
never present at a level after it disappears. We define h{e) and h'{e) to be the recursion levels at 
which the edge e appears and disappears, respectively. 

For every edge e and every recursion level i at which it is present, we let U{e, i) denote the set 
of vertices V of the edge-contracted graph containing its endpoints. If h{e) < i < h'{e), then we let 
W{e,i) denote the set of vertices Vj output by ImpStarDecomposition that contains the endpoints 
of e. 

Recall that p(j) denote the index-mapping of the vertex set Vj in the star decomposition. For 
each index z G {0, 1, . . . , t — 1}, let /j = {j £ {1, 2, . . . , A;} | p{j) = i}. For a vertex subset U V, 
let AS(f7) denote the average stretch that the algorithm guarantees for the edges of E[U). Let 



19 



TS(C/) = AS{U) ■ \E{U)\. Then by Lemma I^TTl the following recursive formula applies. 




where we recall (3 = e = (21og4/3(n + 32))^-"^. 

For every edge e and for every h{e) < i < h'{e), let TTi{e) denote the index- mapping of the 
component W{e, i) in the invocation of Algorithm ImpConeDecomp on recursion level i. For every 
index p G {0, . . . ,t — 1}, define the variable lp{e) as follows 

lp{e) = \{h{e)<i<h'{e)\TTi{e)=p}\ . 

For a fixed edge e and an index p £ {0, . . . , t — 1}, every h{e) < i < h'{e) such that TTi{e) = p 
reflects a contribution of 0{t/ 0) • log*-^^^^/*(m + 1) to the right term in ()12() . Summing p over 
{0, 1, . . . , t — 1}, we obtain 

t-1 

^0(t//3)Zp(e)log(*'+i)/*(m + l). 

p=0 

In a few moments, we will prove that 

t-1 

^5^/p(e)logP/*(m + l) < 0(m log2m), (13) 
e p=0 

which implies that the sum of the contributions of all edges e in levels h{e) < i < h'{e) to the right 
term in ((T^ is 

o(j^-m log^+^/* fhj =o(t-m log^+^Z* 

As vol {Vj) counts the internal edges of Vj as well as its boundary edges, we must also account 
for the contribution of each edge e at level h'{e). At this level, it will be counted twice — once in 
each component containing one of its endpoints. Thus, at this stage, it contributes a factor of at 
most 0{{t/ (3) ■ logm) to the sum TS(y). Therefore all edges e £ E contribute an additional factor 
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of 0{t ■ mlog^ fh). Summing over all the edges, we find that all the contributions to the right term 
in H12|) sum to at most 

o(t-m log^+^/* mj . 

Also, every h{e) < i < h'{e) such that the edge e belongs to the central component Vq of the 
star decomposition, reflects a contribution of O(logm) to the left term in (|12() . Since there are at 
most O(logm) such is, it follows that the contribution of the left term in (|12jl to TS(y) sums up 
to an additive term of 0(log^ fh) for every single edge, and (mlog^ fh) for all edges. 

It follows that TS(V^) = 0{t ■ mlog^"*"^/* m). This is optimized by setting t = log log m, obtain- 
ing the desired upper bound of 0(log^n • log log n) on the average stretch AS(y) guaranteed by 
Algorithm ImpLowStretchTree. 

We now return to the proof of We first note that /o(e) is at most O(logm) for every edge 
e. We then observe that for each index p > and each h{e) < i < h'{e) such that 7rj(e) = p, 
vol(E(C/(e,i))) /vo\{E{W{e,i))) is at least 2^"'^""^ (by Lemma lO ((TT])). For h{e) < i < h'{e), 
let gi{e) = vol{E{U{e,i + 1))) /vol {E{W{e,i))). We then have 

n (2'-^^*^)^^^^^ < - n 9^ie) , 

l<p<t-l h{e)<i<h'{e) 

hence X]p=i ^pI^) log^/* m < log m + Yjh{e)<i<h'{e) log 9i{^)- We will next prove that 

X] X] log5fj(e) < mlogm, (14) 

e h{e)<i<h'{e) 

which implies (fT^ . 

Let Ei denote the set of edges present at recursion level i. For every edge e £ Ei such that 
i < h'{e), we have 

9^{e') = 5i(e)vol {E{W{e, i))) = vol {E{U{e, i + 1))) , 

e'eE{W{e,i)) 

and so X]ee-B,-i</i'(e) ~ vol(i?j+i). As each edge is present in at most O(logm) recursion 
depths, vol (Ei) < fhlogfh, which proves (fTl|) . □ 

6 Conclusion 

At the beginning of the paper, we pointed out that the definition of stretch used in this paper 
differs slightly from that used by Alon, Karp, Peleg and West If one is willing to accept a 
longer running time, then this problem is easily remedied as shown in Subsection 11.11 If one is 
willing to accept a bound of O(log^n) on the stretch, then one can extend our analysis to show 
that the natural randomized variant LowStretchTree, in which one chooses the radii of the balls 
and cones at random, works. 

A natural open question is whether one can improve the stretch bound from O (log^ n log log n) 
to O(logn). Algorithmically, it is also desirable to improve the running time of the algorithm to 
0(m log n). If we can successfully achieve both improvements, then we can use the Spielman-Teng 
solver to solve planar diagonally dominant linear systems in 0(n log nlog(l/e)) time. 
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As the average stretch'^ of any spanning tree in a weighted connected graph is 1^(1), our low- 
stretch tree algorithm also provides an O (log^ n log log n) -approximation to the optimization prob- 
lem of finding the spanning tree with the lowest average stretch. It remains open (a) whether our 
algorithm has a better approximation ratio and (b) whether one can in polynomial time find a 
spanning tree with better approximation ratio, e.g., O(logn) or even 0(1). 
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